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Summary 


The performance specifications of any actuator are quantified in terms of an exhaustive 
list of parameters such as bandwidth, output control authority, etc. Flow-control applications 
benefit from a known actuator frequency response function that relates the input voltage to the 
output property of interest (e.g., maximum velocity, volumetric flow rate, momentum flux, etc.). 
Clearly, the required performance metrics are application specific, and methods are needed to 
achieve the optimal design of these devices. Design and optimization studies have been 
conducted for piezoelectric cantilever-type flow control actuators, but the modeling issues are 
simpler compared to synthetic jets . 1 Here, lumped element modeling (LEM) is combined with 
equivalent circuit representations to estimate the nonlinear dynamic response of a synthetic jet as 
a function of device dimensions, material properties, and external flow conditions. - ’ These 
models provide reasonable agreement between predicted and measured frequency response 
functions and thus are suitable for use as design tools. In this work, we have developed a 
Matlab-based design optimization tool for piezoelectric synthetic jet actuators based on the 
lumped element models mentioned above. Significant improvements were achieved by 
optimizing the piezoceramic diaphragm dimensions. Synthetic-jet actuators were fabricated and 
benchtop tested to fully document their behavior and validate a companion optimization effort. 
It is hoped that the tool developed from this investigation will assist in the design and 
deployment of these actuators. 
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Optimization of Synthetic Jet Actuators 


Problem Description 

Synthetic or “zero-net mass flux” jets 4 are commonly used as flow-control actuators in a 
wide spectrum of applications including jet vectoring , 5 separation control , 6 ' 7 ' 8 and boundary layer 
control . 9 ’ 10 The performance specifications of any actuator are quantified in terms of an 
exhaustive list of parameters such as bandwidth, control authority, etc. Flow-control 
applications require a known actuator frequency response function that relates the input voltage 
to the output property of interest (e.g., maximum velocity, volumetric flow rate, momentum flux, 
etc.). Clearly, the required performance metrics are application specific, and methods are needed 
to achieve the optimal design of these devices. This work summarizes lumped element modeling 
(LEM) of synthetic jet actuator and presents results for the optimization of synthetic jet actuators 
based on LEM, namely the optimization of the cavity and orifice dimensions and the 
piezoelectric composite diaphragm. 

Objectives 

The objectives of this report are: 

■ To summarize the use of lumped element modeling as a design tool for the 
optimization of synthetic jet actuators 

■ To optimize the cavity and orifice dimensions of the actuator for a given piezoelectric 
diaphragm 

■ To optimize unimorph and bimorph piezoelectric diaphragm for various 
configurations 

■ To develop a Matlab-based optimal design synthesis tool 
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Approach 


We have investigated various actuators for active flow-control: oscillating flaps and 
synthetic-jet actuators. In addition, lumped element modeling has been successfully applied to 
piezoelectric-driven synthetic jet actuators as a design tool. In this research, synthetic jet 
actuators are modeled and optimized based on lumped element modeling. 11 ' 12 ’ 13 

Lumped Element Modeling 

A cross-sectional schematic and corresponding equivalent circuit representation for a 
typical piezoelectric-driven synthetic jet are shown in Figure 1, using LEM. The details of the 
circuit parameter estimation techniques, assumptions, and limitations are given by Gallas et al.~ 
The structure of the equivalent circuit is explained as follows. A harmonic voltage V ac is applied 

across the piezoceramic to create an effective acoustic pressure that drives the diaphragm into 
motion. This represents a conversion from electrical energy to acoustic energy and is 
represented by an ideal transfonner possessing a turns ratio <j ) a . The motion of the diaphragm 

(i.e., volume velocity, Q) can store potential energy via compressibility effects in the cavity ( Q c ) 
and/or can store kinetic energy via oscillatory flow through the orifice ( Q out ). Physically, this is 
represented as a volume velocity divider, Q = Q C + Q out . 

There are several simplifying assumptions employed in the current model. First, the 
synthetic jet is assumed to exhaust into a semi-infinite quiescent air medium. In practice, these 
devices interact with a boundary layer that greatly alters the jet-exit velocity profile and thus the 
total orifice impedance. 13 The corresponding differences in the frequency response functions of 
a synthetic jet actuator exhausting into a quiescent medium versus one exhausting into a 
boundary layer means that the bench-top calibration of these devices is insufficient to accurately 
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estimate the volume velocity or momentum flux exhausted into a cross flow for a given 
excitation voltage. Second, compressibility effects in the orifice, but not in the cavity, are 
neglected. The incorporation of these effects into the lumped element model is an ongoing 
research area in the engine nacelle acoustic liner community. 16 We have completed some 
extensions of the lumped element modeling technique to handle a grazing boundary layer. The 
details are beyond the scope of this report, but the interested reader is referred to Gallas et 
al. 17,18,3 



Figure 1: Schematic representation and equivalent circuit model of a piezoelectric-driven synthetic jet. 

The frequency response function of the volume flow rate through the orifice per applied 
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voltage for the equivalent circuit shown in Figure 1 is a four-pole, single-zero dynamic system, 

Qq,A s ) m 

T r ( \ 4, 3, 2, 

r ac [s) a 4 s + a i s + a 2 s + a,. s' + 1 

where 


6 


a i - C aD ( R aO + R aN + R al) ) + C aC ( R aO + R aN ) > 

a i = C aD (M aRad + M aN + M aD ) + C aC ( M aRad + M aN ) + C aC C aD R aD ( R a0 + R aN ) , 

{ 2 } 

a 2 = C aC C aD \_M aD ( R aO + R aN ) + aRad + ^ aN ) R aD ] ’ aI1 ^ 
a 4 = ^aC^-aD^ aD aRad + ^ aN )• 

To first order, the coefficients in Eq. {2} are constants determined via simple algebraic 
expressions as a function of geometry and material properties. This model includes one 
empirical constant, R aD , that represents the structural damping of the piezoelectric composite. In 
general, some of the coefficients exhibit frequency and amplitude dependence (i.e., due to 
nonlinear effects). Specifically, R a0 is a non-linear orifice resistance that is proportional to the 

volume velocity, Q out (s) . 

For a dc voltage (5 = 0), the volume velocity is zero. At low frequencies (5 — » 0), the 
volume velocity is proportional to j<x>d a V ac . At high frequencies (5 — » x ), 

^ ^ out ^ Cl j 1 

7T = 

and the output attenuates at 60 dB/decade. 

The four-pole system in Eq. { 1 } possesses two resonant frequencies, f x and f 2 > , that 

are related to the short-circuit piezoelectric diaphragm natural frequency f D , 


f D = ir- 


2 n\M aD C aD 


and the Helmholtz resonator frequency f H , 


{4} 


by the equality 


fn= — 


2n ]j (M oN + M aRad ) C aC 


{5} 
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{ 6 } 


/i/ 2 = fafn ■ 


From Eq. {6}, it is obvious that the cavity volume and the orifice dimensions, as well as the 
piezoelectric-diaphragm characteristics determine the dynamic response of the synthetic jet. 

Gallas et al. have experimentally validated the lumped element model for two different 
prototypical synthetic jet actuators using phase-locked Laser-Doppler Velocimetry. The 
amplitude of the piezoelectric excitation voltage was 25 V in all cases. The comparison between 
the full nonlinear model prediction and the experiments are shown in Figure 2 and Figure 3. In 
both figures, the centerline velocity magnitude is plotted as a function of frequency. For the 
lumped element model predictions, the centerline velocity was estimated by modeling the flow in 
the orifice as flow in a circular duct driven by an oscillatory pressure gradient. The first case (I, 
Figure 2) clearly illustrates the two resonant frequencies of the coupled oscillator. The second 
case (II, Figure 3) corresponds to a system possessing a single dominant peak. The lower peak at 
f x ~ 3 1 5 Hz is heavily damped in this case due to the frequency dependent nonlinear orifice 

resistance term R a0 . In both cases, there is sufficient agreement between prediction and 

experiment to justify the employment of LEM as a design and optimization tool. It is important 
to note that the underlying assumptions used to derive each lumped element limit the applicable 
frequency range of the corresponding element from dc to some upper limiting frequency (see, for 
example, Rossi 14 ). 
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Figure 2: Comparison between the lumped element model and experiment for Case I. 2 



Figure 3: Comparison between the lumped element model and experiment for Case II. 2 

The lumped element model presented above is thus a powerful design tool that enables 
the multi-energy domain dynamic modeling of a synthetic jet actuator. This model is therefore 
used as a vehicle to enable optimization. The optimization problem is solved via MATLAB’s 
optimization toolbox. 

Cavity and Orifice Optimization 

Details of the modeling procedure are presented in Gallas, and optimization results are 
presented in Gallas et al. 11,13 First, a constrained optimization of the cavity volume and orifice 
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dimensions of the two baseline synthetic jets presented above, each with a given piezoelectric 
diaphragm, has been performed. The goal is to improve the performance of the nominal designs 
presented in Figure 2 (Case I) and Figure 3 (Case II). The objective function, constraints, and 
design variables are summarized in Table 1. Specifically, Case I is optimized using one cost 
function, while Case II is optimized using a different cost/objective function. 


T a ble 1: Summary of optimization problem for a synthetic jet actua tor. 


Objective 

Constraints 

Variables 

= integrated 

centerline 

velocity 

over a 

desired 

frequency 

range 

• Upper/Lower bounds on variables: 
Case I: 

f 0.125 < « 0 < 2 

< 0.5 < Z < 20 
[ll52<F 0 < 47,640 

Case II: 

f 0.125 < a 0 < 2 

< 0.38 < Z < 20 
[780 <V 0 <21,270 

• Orifice aspect ratio L/a 0 > 1 

• Fixed input voltage V ac = 25 V 

a 0 [mm] 
L [mm] 
V 0 [nun 3 ] 

w 0 (/ 2 ) = 

centerline 
velocity at a 
fixed 
frequency 

fi 

• Upper/Lower bounds on variables: 
Case II: 

f 0.125 < a 0 < 2 
j 0.38 < Z < 20 
[780 <V 0 <21,270 

• Natural frequency 
f 2 = 820 Hz 

• Orifice aspect ratio L/a 0 > 1 

• Fixed input voltage V ac = 25 V 

a 0 [mm] 
L [mm] 
V 0 [mm 3 ] 


The first cost function employed maximizes the integrated centerline velocity over the 
entire frequency range, f u 0 (f)df where the upper limit of integration is 3000 Hz and 1500 

J 0 


Hz for Case I and II, respectively. The motivation for such an objective function is to increase 
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the broadband response of the actuator. The constraints on the orifice radius are motivated by 
device manufacturability and flow perturbation concerns, while the minimum orifice length 
constraint is driven by the requirement that the orifice plate be rigid. The volume range is 
dictated by size limitations. In addition, constraints are imposed on the orifice aspect ratio (since 
the current lumped element model is not valid for small L/ a 0 ) and by fixing the piezoelectric 

drive voltage. Figure 4 shows the resulting optimized frequency response function compared to 
the nominal response for Case I. The frequency response function increased over the entire 
frequency range by decreasing V () , a 0 , and L . In addition, the first resonant peak is heavily 
damped in the optimized response. 



Frequency [Hz] 

Figure 4: Optimization of Case I, maximizing the overall output. 

The results using a similar approach for Case II are shown in Figure 5. Again, the 

frequency response function increased over nearly the entire frequency range by decreasing V 0 , 
a 0 , and L . In addition, the flat portion of the response function is increased and the second 
resonant peak is moved to a higher frequency. This actuator design is useful for applications that 
require a flat broadband response. 
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Figure 5: Optimization of Case II, maximizing the overall output. 

The second cost function employed maximizes the centerline velocity at the second 

resonant frequency of the system, u 0 (f 2 ) . From a practical standpoint, this optimization is 

useful for applications requiring high actuation authority over a narrow frequency range. The 
constraints on the orifice and cavity geometry, as well as the aspect ratio and drive voltage are 
similar to those outlined above. A new equality constraint, however, is placed on the location of 
the second resonant frequency. The results of this optimization for Case II are shown in Figure 
6. The second resonant peak is increased as desired (by ~ 35%). In addition, the broadband 
response, especially at the lower frequencies, is increased with respect to the nominal case. 



Figure 6: Optimization of Case II, maximizing the second peak. 
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While the optimization studies in this section indicate the promise of improved 
performance, additional gains can likely be achieved by optimizing the piezoelectric composite 
diaphragm. 

Piezoelectric Diaphragm Optimization 

The piezoelectric diaphragm optimization is described in this section. The configuration 

here is the standard unimorph inner-disc piezoceramic patch bonded to a metal shim. The basic 
modeling concepts of the piezoelectric diaphragm are explained as follows. Ideally, a 
piezoelectric actuator is a linear, conservative, reciprocal transducer . 14 The piezoelectric 
composite deforms in response to both an applied ac voltage and a differential pressure. The 
lumped piezoelectric coupling equations for the transduction model are 19 



C 


aD 


d„ 


d a ]\P 

c EF \\r, 


{ 7 } 


where AV is the volume displaced by the plate due to the application of a differential pressure 
P and/or voltage V ac , q is the charge stored on the piezoelectric electrodes, C EF is the electrical 
free capacitance of the piezoelectric material, C aD is the short-circuit acoustic compliance of the 
plate, and d a is the effective acoustic piezoelectric coefficient. 

The acoustic mass M aD is determined by equating the lumped kinetic energy of the 
vibrating diaphragm expressed in acoustic conjugate power variables to the total kinetic energy. 
The determination of the lumped element parameters for this two-port model requires the 
solution of the transverse static deflection field as a function of pressure and voltage loading. 
For this study, linear laminated plate theory is used to solve for the transverse static deflection, 
and an optimization scheme is then implemented using this two-port model. The details of the 
composite plate model are presented in Prasad et al . 19 The optimization objective function 
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maximizes the achievable displaced volume per applied voltage, d a = ^V\ p=Q /V ac max , compared 

to the nominal values for Cases I and II. The optimization procedure (design variables, 
constraints, cost function, etc.) is described in Gallas et al. 11 . The results are presented in Table 
2. The optimum design slightly improves (by ~6%) the volume displacement achieved for the 
nominal Case I, while a 52% improvement is achieved compared to the nominal Case II. 

Table 2: Optimization results of piezoelectric diaphragms with an inner-disc conf iguration. 


Table 2: Optimization results of piezoelectric diaphragm 
configurations: Case I, =11.5 mm, / D =2114Hz. Max(Q). no 
minimum gauge for t p and t s 


Original 

Optimum 

Unimorph 

Optimum 

Bimorph 

R 3 (mm) 

11.5 

11.5 

11.5 

R 2 (mm) 

Na 

11.5 

11.5 

R t (mm) 

10.00 

11.39 

11.31 

t s (mm)/(mil) 

0.150/5.9 

0.024/0.9 

0.024/0.9 

t p ( mm) /(mil ) 

0.08/3.1 

0.243 / 9.6 

0.135/5.3 

Id (Hz) 

2114 

2114 

2114 

Q xIO 10 (m 3 ) 

52.3 

155.7 

267.2 

V max (V) 

95 

288 

159 
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Table 3: Optimization results of piezoelectric diaphragms with an inner-disc conf iguration. 


Table 3: Optimization results of piezoelectric diaphragm 
configurations: Case II, ^=18.5 mm, / fl =632Hz. Max(Q). no 
minimum gauge for t and t s 



Original 

Optimum 

Unimorph 

Optimum 

Bimorph 

R 3 (mm) 

18.5 

18.5 

18.5 

R 2 (mm) 

na 

18.5 

18.5 

R t (mm) 

12.50 

18.32 

18.32 

t s (mm)/(mil) 

0.10/3.9 

0.018/0.7 

0.0062/0.2 

t p (mm)/(mil) 

0.11/4.3 

0.188/7.4 

0.11/4.3 

fn(Hz) 

632 

632 

632 

QxlO ~ 10 (m 3 ) 

561 

1348 

2460 

V max (V) 

130 

222 

130 


Table 4: Material properties. 


Property 

Piezoceramic (PZT-5A) 

Shim (Brass) 

Elastic Modulus (Pa) 

6.3xl0 10 

8.963xl0 10 

Poisson’s Ratio 

0.33 

0.324 

Density (kg/m 3 ) 

7700 

8700 

Rel. Dielectric Constant 

1750 

- 

d n (m/V) 

-1.75xlO' 10 

- 


In addition, a study to find the optimal piezoelectric composite circular plate 
configuration and design for maximized volume displacement and bandwidth has been 
conducted. Their tradeoff and practical manufacturing requirements are considered via 
optimization techniques. The results for brass shim and PZT-5A (Table 3) combinations are 
summarized below. 
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The analytical composite plate model has been verified for a generic design via a 
commercial FEA code (see Figure 7). 


unimorph 


bimorph 



r(m) 


r (m) 


Design 


Ri (mm) 

R : (mm) 

R 3 (mm) 

t p (mm) 

t s (mm) 

16.89 

17.58 

18.50 

0.123 

0.081 


Figure 7: Validation of analytical model via Abaqus FEM code. The vertical deflection per unit applied 

voltage is shown plotted vs. radius. 


♦♦♦ An optimization study has been conducted to detennine the optimum geometry of 
various unimorph and bimorph configurations (Figure 8). The objective function 
chosen was to maximize the static (i.e., dc) volume displaced by the diaphragm, Q 
(at the maximum operational voltage of the device), subject to constraints on the 
minimum allowed natural frequency ( f D ) of the diaphragm and other geometric 

constraints. The resulting designs tended towards full coverage of the shim by the 
inner piezoceramic layer, that is R U R 2 — > R 3 . 
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(a) Circular piezoceramic unimorph composite with 
oppositely polarized inner disc and outer ring combination. 
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(b) Typical AC voltage application on unimorph. 
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(d) Typical AC voltage application on bimorph. 


(c) Circular piezoceramic bimorph composite with 
oppositely polarized inner disc and outer ring combination. 


Figure 8 : Circular composite piezoceramic plate configurations for synthetic jet optimization. Arrow tips 
represent the polarization direction and denotes the positive upward electric field applied through a patch. 


♦♦♦ Compared to the Case I and II results (Table 2 and Table 3) reported in Gallas et 
al., 11 an optimized unimorph configuration increases the volume displacement 
by factors of 3 and 2.4, respectively, for the same natural frequency and R 3 if no 

practical lower bound is applied for the shim and piezoceramic layer thickness. 
Similarly, an optimized bimorph configuration increases the volume 
displacement by factors of 5.1 and 4.4, respectively. This increase translates 
directly into increased velocity output of a synthetic jet. 

♦♦♦ An optimum bimorph was designed to achieve a target volume displacement 
while simultaneously maximizing f D with estimated lower constraints imposed 

on t s and t due to manufacturing limitations, 3 and 5 mil, respectively. The 
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results, which are directly dependent on these constraints, are summarized in 
Table 5 below and indicate that a powerful piezoelectric synthetic jet with 
reasonable bandwidth can be achieved via optimization. Our model suggests f D 

as much as 1660 Hz. which may be increased further if the thinner shim thickness 
can be manufactured. 

Table 5: Optimization results of bimorph composite with £>=430x10 10 m 3 and max(/ fl ) 


with several minimum gauge for t s and 5 mil minimum gauge of t . 


Design 

t s 

3 mil 

t s 

1 mil 

t s 

0.5 mil 

t s 

0.3 mil 

t s 

0.03 mil 

R 3 (mm) 

14.5 

14.5 

14.5 

14.5 

14.5 

R 2 (mm) 

14.5 

14.5 

14.5 

14.5 

14.5 

R (mm) 

13.60 

14.35 

14.35 

14.48 

14.48 

t s (mm)/(mil) 

0.0762 / 3 

0.0254 / 1 

0.0127/0.5 

0.00762/0.3 

0.002/0.08 

t p (mm)/(mil) 

0.1434/5.6 

0.228/8.9 

0.2359/9.3 

0.2447/9.6 

0.2486/9.8 

fn(Hz) 

1659 

2180 

2228 

2278 

2298 

Q xl O' 10 (m 3 ) 

430 

430 

430 

430 

430 

V max (V) 

169 

269 

279 

289 

294 


Conclusions 

The optimization of piezoelectric-driven synthetic jet actuators based on LEM has been 
carried out. It has been shown that LEM is a viable tool to optimally design these devices for 
candidate applications. To simplify the problem, the current study has split the optimization 
problem into two parts by separately optimizing the (1) cavity volume and orifice dimensions 
and (2) the piezoelectric diaphragm. Significant enhancements in perfonnance are obtained. A 
Matlab-based design optimization tool has been developed. Future work will add a graphical 
user interface to the Matlab code and will develop a corresponding user’s guide. 
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